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ABSTRACT 

A  closed  solutiom^ars-feo^tt  developed  to  calculate  the  critical 
buckling  load  of  a  uniformly  loaded  sphere  which  yields  results  that 
are  in  close  agreement  with  test.  This  procedure  permits  the  cal¬ 
culation  of  buckling  loads  of  spherical  shells  with  varying  shell 
segments,  varying  loading  (provided  axisymmetric  and  radial),  and 
a  varying  thickness.  ,  ^ 
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SYMBOLS 
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€ 
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<f> 
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V 

I 

A,  B 
w 

v 
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radius 

shell  thickness 
modulus  of  elasticity 
Poisson  ratio 
external  load 
strain,  meridional 
strain,  circumferential 
angle,  meridional 
angle,  circumferential 
denotes  a  change 
external  work 
internal  work 
total  energy 
arbitrary  constants 
deflection,  radial 
deflection,  circumferential 
critical  stress  (PR/2h) 
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INTRODUCTION 

Numerous  articles  have  been  written  with  the  intent  of  trying  to 
obtain  a  solution  for  the  buckling  of  spherical  shells.  Some  solutions 
have  been  found  but  are  based  upon  small  angle  approximations  and 
uniform  load. 

The  purpose  of  this  report  is  to  try  to  obtain  a  solution  without  the 
use  of  small  angle  approximations  or  the  necessity  of  uniform  load  and 
be  applicable  to  any  shell  segment. 

DISCUSSION 

The  problem  under  consideration  is  that  of  elastic  buckling  of  a 
thin  spherical  shell  subjected  to  a  uniform  load. 


The  principle  of  virtual  work  shall  be  employed.  Several  authors 
(Ref.  1,  Z,  and  3)  give  the  form  of  the  energy  expressions.  However, 
modifications  are  required  to  be  acceptable  for  use  here. 
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The  total  energy  of  a  system  consists  of  three  parts: 

1.  Energy  due  to  bending  the  membrane 

2.  Energy  due  to  stretching  the  membrane 

3.  External  work 

Only  the  last  of  the  three  energy  equations  shall  differ  from  those 
given  in  References  1,  2,  and  3. 

The  energy  that  is  contributed  by  bending  shall  be  neglected  as  it 
is  but  a  small  percentage  of  the  total  amount  in  the  shell.  However, 
because  of  bending  there  is  a  change  in  surface  area. 

The  stress  distribution  shall  be  assumed  such  that  the  transverse 
cross  sections  which  were  originally  plane  and  normal  to  the  centerline 
remain  plane  after  bending.  Let  ab  and  cd  be  two  adjacent  cross 
sections  of  the  shell  (see  Fig.  2)  and  d(f>  be  the  small  angle  between 
ab  and  cd  before  bending.  As  the  load  is  applied,  bending  of  the 
cross  section  occurs  and  section  cd  rotates  with  respect  to  ab  about 
the  neutral  axis  nn  by  a  small  angle  Ad<£.  The  angle  is  positive  if 
the  initial  curvature  of  the  shell  is  redt  ced.  Because  of  the  rotation 
the  fibers  on  the  outer  surface  are  compressed  and  the  fibers  on  the 
inner  surface  are  put  in  tension.  Let  y  be  the  distance  from  the 
centroidal  axis  perpendicular  to  the  plane  of  bending,  taken  positive 
in  the  direction  toward  the  center  of  curvature  of  the  centerline  of  the 
element  and  e  is  the  distance  to  the  neutral  axis  nn  from  the 
centroidal  axis.  It  can  be  seen  that  the  extension  of  any  fiber  during 
bending  is  (y  —  e)Ad $  and  that  the  unit  elongation  of  the  fiber  is 

(y-e)Adj. 

<t>  <R-y)4 

The  stress  distribution  is  no  longer  linear  but  follow's  a  hyperbolic 
law.  Assuming  that  the  neutral  axis  and  the  centroidal  axis  are  in  the 
same  plane,  that  the  greatest  stress  is  produced  in  the  innermost 
fibers,  and  that  the  thickness  is  small  in  comparison  with  the  radius, 
we  can  write 

hAd<£> 

*4>  =  ~  2 


Note:  The  minus  sign  indicates  inner  surface. 
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Similarly, 


h&d0 

‘0 =  “  zRde 

Thus,  the  decrease  in  the  angles  during  bending  is 
&d<f)  = 


Ad0 


ZtyRd 9 
h 


(2) 

(la) 

(2a) 


The  original  area  over  which  the  external  load  was  applied  is 
equal  to 

dA  =  Rd0ad£ 

where 
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The  internal  work  neglecting  bending  is 


V  =  R2 


sin  <$d(pdd 


Th^  magnitude  of  the  compressive  stresses  in  the  shell  we  shall 

call 


N 

If  we  assume  that  the  quantities  v  and  w  represent  components  of 
small  displacements  during  buckling  from  the  compressed  spherical 
shell  form,  then  the  N<£  and  Ng  compressive  forces  differ  but  little 
from  the  compressive  forces  77,  thus 

+  N'a 


and 


Ne  =  N0  +  N'e 

Since  the  external  load  P  produced  the  N  forces  in  the  shell 
with  the  original  deformation  w,  the  external  work  shall  include 
only  the  small  increase  as  will  the  internal  work. 

Therefore,  the  internal  work  becomes 

v  =  rZ  fpfe  (N^°  +  N'6*ty  sin  0d<pdd  (6) 

However, 

Eh  , 

N'o  s  J^)K°  *  “ 


+  u‘*) 


and 


Substituting  into  equat.on  6 


v  *  +  ***  +  2m<t>*0)  sin  ti>d(^d6  (6a> 
Then  eqitation  5a  becomes 

T  =  +  ~^~P  (**£  +  sin  (5b) 

The  total  energy  of  the  system,  T  +  V,  is  equal  to  a  constant,  I. 

T  +  V  =  I  (7) 

This  total  potential  energy  does  not  change  when  the  structure 
passes  from  an  equilibrium  position  to  an  infinitesimal,  adjacent 
position. 

Therefore,  the  structure  will  be  in  a  stable  position  if  the  total 
potential  energy  I  is  a  stationary  value  and  is  a  minimum. 

According  to  the  mathematical  rule,  I  will  be  a  minimum  if  the 
second  variation  6ZI  is  positive  for  any  virtual  displacement. 

The  Ritz  method  is  very  general  and  applies  to  all  problems  in 
mechanics  and  physics  which  may  be  considered  as  problems  of  the 
calculus  of  variations. 

If  the  deflection  w  is  expressed  in  the  form 

W  =  +  a2^2^2  +  •  -  •  +  *n<i>n6 n 

where  the  total  function  satisfies  the  boundary  conditions.  The  co¬ 
efficients  a  are  arbitrarily  chosen  parameters. 

If  this  deflection  curve  is  introduced  into  the  energy  equations, 
we  arrive  at  an  expression  for  the  total  potential  energy  I  as  a 
function  of  the  n  parameters  a.  If  w  is  to  be  regarded  as  a 
solution  to  the  extremum  problem,  it  must  satisfy  the  condition  of 

extremum 

T  +  V  =  I  =  stationary 

The  problem  is  then  a  maximum  —  minimum  of  calculus. 


By  taking  the  partial  derivative  of  the  total  potential  with  respect 
to  each  of  the  arbitrary  parameters  and  set  equal  to  zero,  we  arrive 
at  a  system  of  n  homogeneous  linear  equations  from  which  the 
parameters  a  are  to  be  determined.  This  system  of  equations  does 
not  have  a  solution  different  from  zero  unless  the  determinant  D 
of  its  coefficient®  is  equal  to  zero.  Therefore, 

D  =  0 


is  the  equation  of  degree  n  in  the  unknown  force  P  and  from  this 
the  stability  condition  from  which  P  can  be  determined.  The  smallest 
root  is  the  P  =  <7. 

Therefore,  the  arbitrary  parameters  a  remain  arbitrary  because 
the  determinant  D  of  their  coefficient  vanishes. 


Success  or  failure  in  applying  the  Ritz  method  to  any  problem  is 
content  upon  the  individuality  of  the  problem  and  the  analysis  may 
become  unnecessarily  lengthy  and  laborious  if  the  series  is  not 
properly  selected. 

The  Ritz  method,  although  giving  values  higher  than  actual,  offers 
the  means  for  the  approximate  solution  to  buckling  problems  in  those 
cases  where  the  exact  solution  may  become  too  difficult  or  is  not 
practicable. 

Timoshenko  provided  a  concept  similar  to  the  Ritz  method. 
However,  whether  Timoshenko  or  the  Ritz  method  is  used,  results 
are  higher  than  actual. 


For  certain  buckling  problems  a  good  approximation  can  be 
obtained  with  only  one  term  in  the  series. 


An  extension  of  the  Ritz  method,  which  was  madar  by  Trefftz  in 
1935,  is  not  included  in  this  report. 


Substituting  the  linear  strain  -  displacement  relationships  (Ref.  4) 


into  equations  6a  and  5b  we  get  the  total  potential  I  of  the?  shell 
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1  -  I ^?fp£  +  -  V-  cot  *) 

4-  v2  cot2  0  4-  Zuv^-^-  cot  sin  0d0d Q 

4-  ~^P~~  f<p  Jg  +  VW  COt  ^  ”  ^w*)  sin  <2>d0d#  (8) 

By  substituting  a  selected  deflection  for  w  and  v  we  can  obtain 
the  following  total  potential 

I  =  XiA2  +  X2AB  +  XjB2  (9) 

where  X  is  a  function  of  (constant,  P,  E,  h,  R,  and  u)  and  A  and  B 
are  arbitrary  constants. 

If  we  take  the  partial  of  equation  9  with  respect  to  A  then  B, 

we  get 


~  =  2X,A  4-  X2B  =  0 

a 

=  X2A  +  ZXj  B  =  0 

<© 


Therefore,  the  determinant  D  is 
D  =  4XjX3  -  X22  =  0 

From  this  we  can  solve  for  the  critical  value  of  the  external  pressure. 

For  a  complete  sphere  with  a  uniform  load  and  u  =  0.3,  we  obtain 
a  critical  stress  as  seen  in  Appendix  A 


a  =  0.  282 


Eh 


which  compares  with  the  test  values  of 


a  =  0.  154 


Eh 

R 


F 


8 


I 


This  calculated  value  is  slightly  higher  than  those  given  by  test 
but  the  procedure  allows  for  the  calculation  of  any  shell  opening, 
loading  (provided  axisymmetric  and  radial)  and  the  use  of  a  varying 
thickness  is  expressed  in  terms  of  the  meridional  angle.  The  select  ed 
form  of  the  deflection  curve  eeries  must  satisfy  the  boundary  conditions 
of  the  specific  problem. 

Appendix  A  is  a  sample  calculation  for  a  two -term  deflection  curve 
for  a  hemisphere. 

Figure  3  shows  the  buckling  curve  for  varying  shell  openings  from 
10  to  90  degrees. 


9 


ww  m  «m  mm-m  mm  www  war  ww  www  h«  ww wwwww;  j 

SaiammSBiiiiKK 

pSsSSsS|i| 

■jM-l'"";--!;:;  “SSS  ihh  KSS 


SH«7«  *■•»«! 

Haawww  •■•I 
■*»*ww  «•*! 

|33Sw5w55Sf 

m«*m 


»:::! 
**•>■ 


SST  - » •2Xa35«3S*#2Sw2SSS] 

Kfc:i::::::::::::::::| 


w#www*  mmmmm  MjaMMM  •  mm'1  «*«««■ 
•  «•■«■  *WW**w  *a*»*  fluid 

.  wwaWwwaMpaawWWwiwi*  WWWW#* 
wwwwwwww  *PwPW#ww«wwwf  •  PfPwww 

mmwwmmimm 


*****  i?*?5*jj*g*  j**|  •  J«*S5«3 


ww  mw  « mmmmm  mmmmm  mmmmm  mmm 
>ww»*a  pww«w*a«*a  wmmwwmwwww  * — 


_ _  -  - _ _ WWWWWWWWWP’ „ 

■  ••««  *MWWW  awwww  wwwww  WWWPP  » 

■  Ml  MMM 

•■*••**••#«<••••  •■«»»  M  . 

wwwww ***** awwww **www  *WW»w  aw>  m. 
mmmmmmrnmmm  mmmmm  rnmmmtt  mmmmm  mm  >wm 
mmmmm  mmmmmmrnmmm  *Wwww**w«w«#  <g<- 


w  »  4  p  ai  ai  M  i  m  mmm  a  www*  an 

••?  fissl  Kg;*  ••••••' 


X 

Ski 
w** 

32 

PwwSppwppiwiajawa 
wapwwwpwwawwpp**  a* 

•  •www*  aw* ww*ww*«  if* 

*  - m  mmmmmmrnmmm  *m» 


a#wwpww< 

•awwwww  aaaaaaai 

mmmmm  mmmmm  mmmmm  mmm 


*:r=«=S::*sEr3=r= 

awwww  mmmm-  wwww*  mm  «| 


»aw«to*»wwa  »-* 

W»*W#WWWWW  •  • 
WWWWW 

WMMiMHlMMMM' 


•  MWM 
n.,  mmmmm  ww 
wwwwww**  aw 
wwwww  1“ 
wwgwtpwww 


•  wa  *w**w**mw*«m«Pw  »•*■«••*•<  wwwwww  wwwww *****  **********  awwww  *****  ***** •»***  *■»*•••' 

••§8? •SfisiiSi!  5???¥5«lS9f?88f*  5ICS8MSM  22221 52S22  22  jSSJ]?1 

•3  Si 

an 

E 

::s 
s 

sf 

1225#  55225 .. 

::::: ::::::::::  :::::::m: :::  :::H : 

mmmmm  mmmmmmrnmmm  wwwww*****  mmm  »W£wwmw ■l!’***  »•*«■#  ** 

r::::::::::::::;;::::::::;;: 

•  Mt»  «®.awP  «••«*  ir»n«  mmmmm  mmm.  tm  *4«mw  mmmmm  mmmmmmrnmmm  mmmwwmwwm*  •  MMWWwa»www  mmmmm  mmmmm  mm< 

**•  aw  ww«ai**Pai*aNi  wwww  *#••*•••#•  wwwww  wwwww  a *•*#* m  mm m ww  •••  ••  iw 
tr-*T  ■■  aw*  a w awwww  *  if*  wwwww**  w«#w*  # www wwwwww  awwwww****  mmmmm  mmmmm  ‘  * 

mmmmm  *4**»***#w ***w«a*#*w ***  awwwww*  wwwwwpppwa  »*«•*•«»*•  mmmmm  wPww*  awwwww****  »**•••#***  .... 

wwwwwapaiN*  aaik  ••*«*«•  •■••••••«*««««**ik««»«««««ai«Mai*  ««*ai*  •••«*••*•*  «rj 

■  iim(  *•••*•*«■•  «*i  aw*  waiw*  - - 

•  •*(<  MaiaiMiiioiiitaMVii  •  m  mmmmm  aw  ««•■««  a>««  «a»«iw - - - 

- aFaiwwwwwa^awarftwiiiM  wwwwwwwwwa  wwaww 

-  «*««b«f«Hiar»ai««fiai*  aiataibai . 

««•««*«  a«*w"  w« w»«  wwwwwjMkai 

i*«««aa  wa*M««aa*  mmmmm  mmmmm  mmmmm  mmmmm  mmmmm  mmmmm  mmmmmmrnmmm  mm-m  •mmmmm  mmmmm  mmmmm 

_ _ _ _ _ _ ^  J#iam»*.wa»»ib«i««i»"*«ai#»«»»apb«»«#»if»««i»aia»ai*ii**«aian»iMia>i»Biai*i‘  -  .•••amwiaiwwwwi  -*• 

«*aa^  *«aawalaw«*  *w>**W  *-«*•»«*  nf*««wa*.*bar*w* - -  - - a - —  - -  - 

- - iwwww  wmwm« •»«•»«•  nawwaiw w«i 

iwwwwawwaa waiaiaaw aw  <** iwarw* vawaiwai 

lawMkP  ftiaipf  ««6m«  w«««r  - 

ssisss::: 

liiHiilll 


■  wwarww  w«ataw  www  ww  waawww  wSSw*  t»w*iw#  wwwww  aaiaiaw  wii  wwa  aw  war  ft 

■  wwwaw*<rwaw»aai**««wawwaiwaw  wPwWwwwwww  iww  '•••••••waaaai 

aintarwa  wwawwawawwawwawwaiww* «•••*•••«*  *»  ••  awwwa ■**■*»•» 
•jMkwaiwawwaawMWMw  aawwwwaiaiiaw  aaaaaaaa*a  wmw  awawa 

kiaaai •****•*«*• *****  »«aaw mmmmmmrnmmm mmm « •••aa* aaaaa »»»*« 

_ _ _  a# www#  «««ai«  wwwww  a«ww»wawww  aawww wai‘  -  .*w#^ma*w»w#i 

mmmmmm  ••  w w-wwwww#  waiwpw  wpaj#w #•!#«#* waiw pw  jaaawiaaai  wa»p.  iMwaa  w*«w>  m».« 
wwwwwi  widw-Pw##**#  wwPww#i#«tiw«wPw«*aww#aBwaw*w«a*  awaiwwarw w«a  •*--  mmmm* 
wawtarS  wwwwa w#aiPa» ##•*•  wwb#w  waawwwwwaaw aaaaaaatka  #wp«w*aiww#  aa«»iaaai 


'•aaaaaa 
.  iaaai*a»-iar  - 
•v  .  WWWPWPW  .*•• 

I  awpppHpwpa  ««<*•  •*! 

-  -*aaa«a«aM>  wan 

—  - - ^  mmmmmm mmm  mmm 

•awwwww# aPwwwwwwwa aw 


::::: 
a*aaa 
«aiwwaww 
.  wa»wwwpp 

'#««#  ***ww*w*S5!!w*3#2 

aaaa aaaaaaaaaa awwwwwr 
•  a««  mmmmmmmrn ■*«  aaaaa 
*•*««  wjaww w  aiw«*iW  pwwwai 
■  ■www  awaiw-wwwwww  #«wjaw 

<*«  WWW  *waiwa««  WWW  iwwww 
««««■  mmmmmmmmmii  mmm wai 


■  *w •* aw  .  ^ 
JWPPW  wajpp 


CSSSSt ; 

«■«•■■•  I 


•«««..  _ 

_ 

aaaaaaaaaaa  ■ 


- aaaa 

. . .  *pp«w  «■«#■••*■■* 

IMP •■■■* Wpwww P»WW* WWW ww 
• — ““•••■■••pp8p* — . 


laaaaaaa 

at::::: 

wwjMMf  aaai 
«#•••■««! 
#  www# 

i 

mmmmm Am w 


mm  mmmmm  aaaaaai 
-  „ ,p*p««pi »pWP 
awwww aaaaaipaa 
■“•^ppppppp  i«r 

!■•  «4*l 

««aa«w*a 
aaaaaa**' 

- PWWWWW# 

wataaaa 

■«■  mmmmj$ 


■^HBiliiH^lisssssssssnssssssnssssassssssssi 

H:«:ss:slssss!«:s::s:w:mcss: 

yMia*«**aaaaa.»a**a  *  m  www  wai^w*  *  im  ««• 


mmmmm  #m< 

— '•■«  ■■ 

•  PPPP  #■#••■■•■#  wbwpwwwwww 
mmmmm mmmmm  mmmmm  mmmmm  mmmmm 
m * a*wa»aaa»aaaa«»aaai 
ip  •■•aaaaaaaaaaa* aaaai 
ip  mmmf+rnmmmm  mmmmmmrnmmm 
Ip  W«'  w»  awapp  awwwawppww 


aaaaa iaa*a 
-#a*waaaaa 
■  ««•  awwpppi 
www* «awaa~ 


.aaaa* aaaaaiaaaa 

8  MW WWW  WWW WWW WWW# 
paaaa •#■«••■#•• 
wmPww.w  aaaaaaaaaa 

IWW  WWWWWWWPWW  ■*•*•  WWWWW 


.  .  aaaaa*  •*•»■***** •••«****«•  w«awaiw«jwl*aiw  taaaaaaaa* 
*  wwwwp  #ww«-w#ppwp  wppp#wa*ww«i*w«pw  *wwpwpwww*  aaaaaaaaa* 
a  «wwwa*waw»*paBPwwww»«w*Wf  wwWw*  ww*  awapawwwww 

2S£2JJ225£2SS;2J5!25S55S25!??5?5!  IZrZTVMZzhiUtZZl 


S.._ .  pppppp# _ 

wwaiww  aaaataaaaa 
wpaHaw  wwwww  wp  www 

- *  wayaww  w««Pwpa»a*ww 

*aaa*»«»Mattfaa«* 

*  wwwww  wwaiww  Vww  w* 
a  wwwww*  ww«#  ana  ww* 


:::l 


wwwww  * www aa*4 ww aaaaaaaaaa 
wwwww  awwwww aM* w w  wwwwMwwwp  - ... 

*w«w* * *aa»aiMh»i aaaawaa*** wwwww wwwww 
•  WWW*  awwwwwajwww  p ^ — - - 


WWWWW  W WWW  * WWWWW’ 


ft 

u 


watSaawi 

„ ^ .aaaa ww **a  w*a#| 

•Ufa  wwwww  wwwww  wwpwwwwi““' 
;%www#  wwwpw wwwww  wwww Www 
•  - w  wwwww  WWWWWWWI 


Mwaaaaai 


a»w  wwa»ww waiaaiaw  wwwww  — 
pw PPPWWWPWJW PWWWWPP WPP 

- a##*!,.,.... 

w  wwwpw a 


P.WWW#  « 
aaaaa  ~ 
•  WPP# 


wwwww  aaawaaaaaf a 
Wwwai*  >aaM««aa«  m 
mmmmm  waiwwwwplawp  aaaaawaaa* 
WWWWW  WWWWW wwwww  wwwww wwwww 
wwwww  aaaaaaaaaa  wwwww  wwwww  .. 
wpwww  *aa waaaaaa aaaaaaaaaaa 
wwwww • wwww wwaaw  aaaaa aaaaa * 
•  a*** wmwwwwwwww  wwwwwwwwwa 
wpwww  wwawppwwwp  pwwwwpwwwp 
wwwww aaaaaaaaaa aaaaaaaaaa  . 
wwwww  aaaaaaaaaa wWwwwpwwwwp 


MMMMM  «••*••**.*» 

_ _ _ — r_.c.aafwii  aaaaaaaaaa 

..'wwww wwwww  wwwwwwwwww p#w wwa aaaa  wwwwwwwwwa 
••■a* aaaaa wwww w  wpwww wpwww* www* wwwww  wwwww 
—  —  wwwwwpwwww  wwmwwmWwWn ■•*a**tfw*« 

waakawaa* aaaaaaaaaa 
-wwwwwwp*  wwa*ww  wwww# 
-awwwww* aaaaaaaaaa 

{■aaaaa wwwww  wwww* 

•aapaa **»••■**«• 
ww.wiww  wa>w«w  mw  www  *  wwwa 
-  pwpww  WWW WWW WWWW 

WWWWW WWW WWW WWWW 

Jwwwwat  aaaawaata* 
'mmmmm  mmmmmmrnmmm 

Bwwwaiaiwww  awaaaaaaaa 
wwwwwwww  aaaaaaaaaa 
wwpwwwwwwwwwwwwwwpwwww wwwww wwww# wwww* 

:::::::::: 

-  — — —  — - »  •**»»  *  WWWWWWWW* 

a wwPww wwwww wwww www www 
-  — ‘'-“‘aaaaaaaaaaa 

- wwww* 

wwwww 
wwwww 
wwwww 
....  .  wwwww 
wwwwwwwaw* 


10 


APPENDIX  A 
SAMPLE  CALCULATIONS 
HEMISPHERE 
P  =  90° 


+  v2cot2©  +  2uv-^  cot®  1  sin®  dodo 

d©  -l 


Assume : 

w  =  A,cososine  +  A£o$3©sin36 

v  =  B,sinosine  +Bsin  3©sin3© 

2 

4^  =  B.cososine  +Bcos3©sin3e  -3 
d®  1  2 

Substituting  and  noting 

/2tt 

/sin2(i)ede=  tt 
Jo 

2n 

sinesin3ede  =  0 

Performing  the  intergration ,  minimize  the 
total  energy  with  respect  to  A^A^B,  ft  Bg 
and  set  in  determinant  form,  we  get 


n 


APPENDIX  A  (Concluded) 
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